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Abstract 

We compute fundamental solutions of homogeneous elliptic differential operators, 
with constant coefficients, on M. n by mean of analytic continuation of distributions. 
The result obtained is valid in any dimension, for any degree and can be extended 
to pseudo differential operators of the same type. 
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Let be P := P(D X ) a pseudo-differential operator, with constant coefficients, 
obtained by mathematical quantization of the function p, i.e : 

pfw = J e^Ko/m- (i) 

Here, and in the following, the notation : 

/(0 = / e-^f(y)dy, 

designs the Fourier transform of /. We note iS(IR n ) the Schwartz space and 
«S'(M n ) the distributions on <S(R n ). A fundamental solution for P is a distribu- 
tion & G <S'(M n ) such that P& = 5. Fundamental solutions play a major role 
in the theory of PDE. For a large overview on this subject, and applications, 
we refer to [1,2]. Apart in some trivial cases, there is few explicit characteri- 
zations of fundamental solutions and most results concern the existence. The 
case of order 3 homogeneous operators, in dimension 3, was treated in [3]. 
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Always for n — 3, the case of certain elliptic homogeneous operators of de- 
gree 4 was also solved in [4]. We are here interested in the case of a definite 
homogenous polynomial pk on R ra , i.e. Pk(£) = £ = and : 

Pk(\0 = \x\ k P k(0- 

Note that k has to be even but we do not impose the spherical symmetry of pk- 
Strictly speaking, it is not necessary to assume that k G N and we can consider 
operators with a conical singularity at the origin. The main motivation is that 
such an operator can generalize the Laplacian (positivity and ellipticity) but 
these operators also play a role in physical optics. To state the main result, 
we introduce the spherical average of g G 5(lR n ) w.r.t. the symbol p : 

A(g)(r)=r n - 1 J g{r9)p k {9)^ d9 . (2) 

Here p k {9) designs the restriction of pk on the sphere, which defines a strictly 
positive function. Clearly, we have A(g) G «S(M + ) and we can extend A(g) by 
for r < to obtain a L 2 function on the line. The main result is : 

Theorem 1 If k < n a fundamental solution & G iS'(M ra ) for P is given by : 
(©, /) = T^y, J p k (6r l fXr6)r n - l - k drd6. 

But, when k > n, we have : 

(6, /) = -C7 fc>B ^^(0) + D ktn J lo g (u)^01(u)du, 

where Ck >n and D k n are universal constants given by Eqs.(8,9). 

The reader can observe the analogy with the Laplacian, see e.g. [1]. In par- 
ticular, one has to distinguish the case of an integrable (resp. non-integrable) 
singularity for Pk{C)~ l for the n-dimensional Lebesgue measure. 

Proof of Theorem 1. Since p k > 0, we define a family p(z) of distributions : 
V/ G S(W) : (P(z),f) = T^J Pk(0 Z f(0dC 

M." 

The r.h.s. is holomorphic for > —1/k and, by continuity, we obtain : 

^7^Jp^rfWt; = f(0) = (5J). (3) 
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The Laurent development of p in z = — 1 can be written : 

— d co 

P(z-1)= + Mo + E ( 4 ) 

j"=-i j"=i 

This point is justified in Lemma 2 below. But, according to Eq.(3), we have : 

]im(p(z-l),P(D)f) = (6,f). 

z — >U 

It is then easy to check that /j,q is a fundamental solution for P. Also, note 
that Eqs.(3,4) provide the set of non-trivial relations : 

P{D)nj = 0, Vj < 0, (5) 

in the sense of distributions of S'(M. n ). The existence of such non-zero //_,-, for 
j < 0, implies the non-uniqueness of fundamental solutions in 5'(lR n ). 

Lemma 2 The distributions p(z — l) are meromorphic on C with poles located 
at rational points Zj t k = — § , j G N. 

Proof. Let g — f ■ Using standard polar coordinates we obtain : 

(p( z - 1), /) = j±- J (r^y-'girey^drde. 

Here Pk(0) is the restriction of pt to S™" 1 . Next, if we define : 

y = (yi,...,y n ) = {rp k {6)^,6), 
we obtain a very elementary formulation : 

{p{z-l)J) = -±- I y^GiyJdy,. 

This new amplitude G is obtained by pullback and integration : 

G(yi) = j y*(g(r,6)r n - 1 \Jy\)dy 2 ...dy n . (6) 
We have G E S(R + ) and G(y l ) = 0(y^- r ) in Vl = 0. Starting from : 



Qk k— 1 

dy\ 



y? = U(kz - fiyt 

3=0 



k(z-l) 



after integrations by parts, we accordingly obtain that : 

<■><* " i). /> = t^t 1 t^-, j rf^GM*, (7) 
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The integral of the r.h.s. is holomorphic in the strip ^ft(z) > — \. Finally, we 
can iterate the previous construction to any order to obtain the result. ■ 

Since z = is a simple pole, the constant term of the Laurent series is : 

For the calculations, we distinguish integrable and non-integrable singularities. 
Case of k < n. The derivative of the rational function provides : 

where 7 is Euler's constant and ip the poly-gamma function of order 0, i.e. : 

7 = lim - — log(m), 

j=l «< 

1>(z) = d z (\o g (r(z)) 1 " U) 



Accordingly, since G vanishes up to the order n — 1 at the origin, we have : 
C k , n J d k yi G( yi )d yi = -C k , n d k n l G(0) = 0. 

On the other side, by derivation of the integral, we find the term : 



log( yi )d k v G( yi )d yi 



(k - 1)1 J 

Since k < n, we can integrate by parts to finally obtain : 



Here, we have replaced the expression for our amplitude, via inversion of our 
diffeomorphism. This proves the first statement of Theorem 1. ■ 

Case of k > n. The term attached to the derivative of the rational function 
is non-zero and provides : 

-C k , n d^G(0) + JL (-l)*^ / lo g ( yi )d k yi G( yi ) dyi . 
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We compute first the derivative of G. Contrary to the case k < n we cannot 
take the limit directly but we will reach the result with the Schwartz kernel 
technic. By Fourier inversion formula we have : 

9^0(0) = ±J e-^{ti) k - l G{u)dudi. 

We can extend the integral w.r.t. u on the whole line by extending G by zero 
for u < 0. Going back to initial coordinates provides : 

The r.h.s. is exactly the derivative of order k — 1 of A(g) defined in Eq.(2). By 
the same technic, the logarithmic integral gives the non-local contribution : 



1 r . ,d k A(g) / NJ 



(k-i) 

In particular the second universal constant of Theorem 1 is : 

1 

Dk ' n = (2tt)" (Jfe-1)!" (9) 

This proves the second assertion of Theorem 1 . ■ 

The main result holds also holds for elliptic pseudo-differential operators with 
homogeneous symbol of degree a e]0, oo[. Using k — E(a) + 1 integrations by 
parts in Eq.(7), where E(x) is the integer part of x, the construction remains 
the same and all formulas can be analytically continued w.r.t. the degree 
a > 0. This continuation is trivial since the functions of Eqs.(8,9), defining 
respectively C n ^ and D n kl are analytic w.r.t. the variable k > 0. 
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